Scaling behavior of an artificial traffic model on scale-free networks 
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In this article, we investigate an artificial traffic model on scale-free networks. Instead of using 
the routing strategy of the shortest path, a generalized routing algorithm is introduced to improve 
the transportation throughput, which is measured by the value of the critical point disjoining the 
free-flow phase and the congested phase. By using the detrended fluctuation analysis, we found 
that the traffic rate fluctuation near the critical point exhibits the l//-type scaling in the power 
spectrum. The simulation results agree very well with the empirical data, thus the present model 
may contribute to the understanding of the underlying mechanism of network traffics. 

PACS numbers: 89.75.Hc, 89.20Hh, 05.10.-a, 89.40.-a 



I. INTRODUCTION 

Complex networks can be used to describe a wide range 
of systems from nature to society, thus there has been a 
quickly growing interest in this area since the discoveries 
on small-world phenomena and scale-free properties 
Q. Due to the increasing importance of large commu- 
nication networks upon which our society survives, the 
dynamical processes taking place upon the underlying 
structures such as traffics of information flow have draw 
more and more attentions from the physical and engi- 
neering communities. Previous studies on understanding 
and controlling traffic congestion on network have a basic 
assumption that the network has a homogeneous struc- 
ture |3|, 1^, [5J, |fj . However, many real communication net- 
works, including the Internet Q and WWW @, display 
the scale-free property. Therefore it is of great interest 
to explore how this highly heterogenous topology affects 
the traffic dynamics, which have brought a number of 
research works in recent years (see the review paper [s| 
and the references therein). 

To improve the transportation efficiency on complex 
networks, Guimera et al. presented a formalism that 
can deal simultaneously with the searching and traffic 
dynamics in parallel transportation systems |lfj|. This 
formalism can be used to optimize network structure un- 
der a local search algorithm through the knowledge of the 
global information of whole network. By using a global 
and dynamical searching algorithm aimed at the short- 
est paths, Zhao et al. provided the theoretical estimates 
of the communication capacity |ll|. Holme proposed a 
much intelligent routing protocol in which the packet can 
detour at obstacle thus guarantee much better perfor- 
mance than the traditional routing strategy which just 
waits at obstacle Since the global information is 

usually unavailable in large-scale networks, Tadic et al. 
proposed a traffic model on a mimic WWW network 0] 
wherein only local information, is available. This work 
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FIG. 1: The order parameter H vs R for BA networks with 
average degree (k) = 4 and of size N = 1500. The estimative 
values of critical point R c of phase transition are 10, 35, 56 
and 35 for f3 =0.3, 0.5, 1 and 1.5, respectively. 



highlighted the relationship of global statistical proper- 
ties and microscopic density fluctuations |14l [15J . In ad- 
dition, Yin et al. [16|)[l3 investigated how the local rout- 
ing protocol affects the traffic condition. The routing 
strategies for the Internet 0] and disordered networks 
[T^ | are also studied. 

Recent empirical studies on communication networks 
have found pervasive evidence of some surprising scal- 
ing properties. One example of such discoveries is that 
the traffic rates of a given link in the Internet (i.e. the 
number of packets or bytes traversing a given link per 
time unit) exhibit the self-similar (or fractal-like) scaling, 
and the multifractal scaling is also found over small time 
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scales |2jJ, El H2, . These empirical studies describe 
pertinent statistical characteristics of temporal dynamics 
of measured traffic rate process and provide ample evi- 
dence that these traces are consistent with long-range 
correlated behavior. Faloutsos et al. concerned the au- 
tonomous system (AS) of the Internet and pointed out for 
the first time that the degree distribution of snapshots in 
the AS level follows a heavy-tailed function that can be 
reasonably approximated by power-law forms |24j. Fol- 
lowing these pioneering works, Percacci et al. found that 
the round trip time (RTT) or Internet delay time can also 
be characterized by a slow power-law decay j2^. This 
finding implies that network inter-connectivity topology 
structure of Internet and transportation delay both have 
a highly heterogenous characteristic. Thus, the assump- 
tion of Poisson process, which had been a major traffic 
model with homogeneous network in the traditional traf- 
fic theory, has clearly lost its validity. This is also the 
very reason why in this article we use scale-free networks 
for simulations. Furthermore, the observation of a phase 
transition between the free-flow phase and the congested 
phase in the Internet traffic is demonstrated by Takayasu 
et al. through both the RTT experiment Ha |53 and 
packet flow fluctuation analysis j^, [2^. However, the 
1/ / -ty p e scaling is only detected near the critical state 
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In this article, we investigate a traffic model on scale- 
free networks, in which packets are routed according to 
the global topological information with a single tunable 
parameter (3. Instead of using the routing strategy based 
on the shortest paths, we give a generalized routing algo- 
rithm to find the so-called efficient path to improve the 
transportation efficiency that is measured by the value 
of the critical point disjoining the free-flow and the con- 
gested phases. Furthermore, by using the detrended fluc- 
tuation analysis, we investigate the statistic properties of 
traffic rate process (i.e. depicted by the number of pack- 
ets per time unit) based on this traffic model, which is in 
good accordance with the empirical data. 
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FIG. 2: (Color online) The critical R c vs (3 for BA network 
with size N = 1225 (up panel) and N = 1500 (down panel). 
Both the simulation (black squares) and analysis (red circles) 
demonstrate that the maximal R c corresponds to j3 ~ 1.0. 
The results are the average over 10 independent runs. 



II. MODEL 

In this paper, we treat all the nodes as both hosts and 
routers [lOl l30|. The model is described as follows: at 
each time step, there are R packets generated in the sys- 
tem, with randomly chosen sources and destinations. It 
is assumed that all the routers have the same capabili- 
ties in delivering and handling information packets, that 
is, at each time step all the nodes can deliver at most 
C packets one step toward their destinations according 
to the routing strategy. We set C = 1 for simplicity. 
A packet, once reaching its destination, is removed from 
the system. We are most interested in the critical value 
R c where a phase transition takes place from free flow 
to congested traffic. This critical value can best reflect 
the maximum capability of a system handling its traffic. 
In particular, for R < i? c , the numbers of created and 



delivered packets are balanced, leading to a steady free 
traffic flow. For R > R c , traffic congestion occurs as the 
number of accumulated packets increases with time, sim- 
ply for that the capacities of nodes for delivering packets 
are limited. To characterize the phase transition, we use 
the following order parameter 



ff( *)=lim^m. 

V ' t^oo R At 



(1) 



where AW = W[t + At) - W(t), with (• ■ • } indicating 
average over time windows of width At, and W(t) is the 
total number of packets in the network at time t. Clearly, 
H equals zero in the free-flow state, and becomes positive 
when R exceeds R c . 

Consider the routing under global protocol [3lJ , where 
the whole structural information is available, and the 
fixed routing algorithm is the most widely used one 
for its obvious advantages in economical and technical 
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FIG. 3: (Color online) The traffic rate process for free (red), 
critical (blue) and congested (black) states with different f3. 
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FIG. 4: (Color online) The detrended fluctuation analysis of 
the traffic rate processes generated by the present model. All 
the data are obtained from the critical state, and the different 
curves represent the cases of different f3 from to 1.9 at step 
0.1. 



costs Actually, the path with shortest length 

is not necessarily the quickest way, considering the pres- 
ence of possible traffic congestion and waiting time along 
the shortest path. Obviously, nodes with larger de- 
gree are more likely to bear traffic congestion, thus a 
packet will by average spend more waiting time to pass 
through a high-degree node. All too often, bypassing 



those high-degree nodes, a packet may reach its desti- 
nation quicker than taking the shortest path. In order 
to find the optimal routing strategy, we define the effi- 
cient path 31]. For any path between nodes i and j as 
P(i -> j) : = i = xq, Xi, ■ ■ -x n -i,x n = j, denote 

n-1 

L(P( J ^j):/?) = E fc ^)' 3 ' ( 2 ) 

i=0 

where k(xi) denotes the degree of the node Xi, and f3 
is a tunable parameter. The efficient path between i 
and j is corresponding to the route that makes the sum 
L(P(i — > j) : (3) minimum. Obviously, L m i n (j3 = 0) 
recovers the traditionally shortest path length. As for 
any pair of source and destination, there may be several 
efficient paths between them. We randomly choose one 
of them and put it into the fixed routing table which is 
followed by all the information packets. 



III. SIMULATION AND ANALYSIS 

Some simulation and analytic results are presented in 
this section. Without particular statement, all the simu- 
lations are based on the Barabasi- Albert (BA) networks 
with average degree (k) = 4. 



A. Phase transition 

The packets handling and delivering capacity of the 
whole network is an index of performance in traffic sys- 
tem, which can be measured by the critical value R c . At 
the critical value R c , a continuous transition will be ob- 
servation from the free-flow phase to the congested phase. 
Fig. 1 shows the order parameter H versus R with dif- 
ferent j3, in which the obvious phase transitions occur. 
For different /3, it is easy to find that the capacities of 
systems are much different. 



B. Optimal routing strategy 

We select packet routing strategy through a tunable 
parameter /?, and the capacity of whole network are much 
different. A natural question arises: Which value of (3 
will lead to the maximal capacity of the network. In 
Fig. 2, we report the simulation results for the critical 
value R c as a function of (3 on BA networks with the size 
N = 1225 and N = 1500, which demonstrate that the 
optimal routing strategy is corresponding to f3 = 1 and 
the size of BA network doesn't affect the optimal value. 
In comparison with the traditional routing strategy (i.e. 
(3 = 0), the capacity R c of the whole network is greatly 
improved more than 10 times without any increase in 
algorithmic complexity. 
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FIG. 5: (Color online) The detrended fluctuation analysis 
of the traffic rate processes generated by the present model. 
All the data are obtained from the free-flow state, and the 
different curves represent the cases of different /3 from to 
1.9 at step 0.1. 



B y e xtending the concept of betweenness centrality 
|34l l35j| to efficient betweenness centrality j3|J , the an- 
alytic results can be obtained according to the Little's 
law (One can find the concept of efficient betweenness 
and the details about the analytical approach in Refs. 
[HEl E3) The analytical results are also shown in 
Fig. 2, which agree very well with the simulations. 



IV. DETRENDED FLUCTUATION ANALYSIS 
ON THE SCALING BEHAVIOR OF THE 
TRAFFIC RATE PROCESS 



critical and congested states, respectively. As time goes 
on, W(t) in the congested state will become too large to 
be plotted together with those in free and critical state. 
Therefore, each curve representing the congested state 
in Fig. 3 is anteriorly divided by a very large number. 
Those numbers are not the same for different simulations, 
and do not have any physical meanings, they are just used 
to make the figure more clear. The behaviors of W(t) 
in the free and congested states are very simple: In the 
former case, it fluctuates slightly around a very low value, 
while in the latter case, it increases linearly. However, the 
time series at the critical point is very complicated, we 
next detect its scaling property by using the detrended 
fluctuation analysis. 



B. Detrended fluctuation analysis 

The empirical study on the self-similar scaling behavior 
of traffic rate process can be firstly found in pioneer ref- 
erences (37, 38] . The autocorrelation function and power 
spectrum are widely used to analyse the the self-similar 
scaling behavior of data [33, H3, EH E2- However, it is 
shown that all the above methods don't work very well 
for the effect of non-stationary, and are less accurate than 
the detrended fluctuation analysis (DFA) O ES EH, 
which has now been accepted as an important time series 
analysis approach and widely used especially for financial 
and biological data [ir l l47 l I4 H li l |5C | . 

The DFA method is based on the idea that a corre- 
lated time series can be mapped to a self-similar process 
by an integration. Therefore, measuring the self-similar 
feature can indirectly tell us the information about the 
correlation properties [43L 0] . Briefly, the description of 
the DFA algorithm involves the following steps. 

(1) Consider a time series Xi, i = 1, • • • ,N, where N is 
the length of this series. Determine the profile 



A. The critical phenomena of the traffic rate 
process 

The empirical studies on communication networks con- 
firm the existence of critical phenomena in traffic rate 
process [2^, |2^. The Congestion duration length dis- 
tribution obeys a power-law form with an approximate 
exponent -1 in critical state, which implies the self- 
similarity scaling exponent (Hurst exponent) H » 1. 
Actually, H < 1/2 stands for an anti-persistent long- 
range correlated process in which a step in one direction 
is preferentially followed by a reversal of direction, while 
H > 1/2 is interpreted as a persistent long-range corre- 
lated process in which a step in one direction is preferen- 
tially followed by another step in the same direction. A 
value of H = 1/2 is interpreted as the ordinary diffusion 
(random walk), in which each step is independent of its 
preceding one. 

In Fig. 3, we report the average number of packets over 
all the nodes, W(t) = W(t)/N, as a time series in free, 



y{i) = ™ fa)]'* = 1 >" ' > N > 



k=l 



where 



1 N 



(3) 



(4) 



(2) Divide profile y(i) into non-overlapping boxes of 
equal size t. 

(3) Calculate the local trend y±i t in each box of size 
t by a least-square fit of the series, and the detrended 
fluctuation function is given as 



Y(k)=y(k)-y lit (k)- 



(5) 



(4) For a given box size t, we calculate the root-mean- 
square fluctuation 



F(t) = 



\ 



1 N 



(6) 
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and repeat the above computation for different box sizes 
t (i.e. different scales) to provide a relationship between 
F and t. If the curve F(t) in a log-log plot displays 
a straight line, then F(t) obeys the power-law form t H 
with H the slope of this line. 

As shown in Fig. 4, the mimic traffic rate process also 
exhibits the self-similar scaling behaviors at the critical 
point of phase transition. The scaling exponents calcu- 
lated with DFA for different are approximate H w 1, 
and the value of /3 has almost no effect on H . This value 
of H implies the 1 / /-type scaling in the power spectrum 
and the long-range correlated behavior in a wide range of 
scales. A very recent empirical study on the traffic rate 
process of a University Ethernet has demonstrated that 
the real Ethernet traffic displays a self-similarity behav- 
ior with scaling exponent « 0.98 [5]|. which agrees well 
with the present result H « f . 

In order to confirm that the critical behavior only ex- 
ists in the critical state (as is shown in some previuos 
empirical studies) (2(J 153, 13 |2i| , we analysis the traffic 
rate processes in the free-flow state. Fig. 5 suggests that 
traffic rate processes in free state are completely different 
from those in the critical state, which exhibit a very weak 
long-range correlated behavior with exponent H = 0.55. 
Clearly, if the transporting time of each packet is ex- 
actly the same as the path length (i.e. a packet does not 
need to wait at any nodes in the path from its source 
to destination), there should be no long-range correla- 
tion and the Hurst exponent will be H = 0.5. In the 
present model, even in the free-flow state, a packet may 
wait some time steps at some nodes before reaching to 
its destination. That is the reason leading to the weak 



long-range correlation. However, the waiting effect in the 
free-flow state is trifling for the Hurst exponent is close 
to 0.5. 



V. CONCLUSION 

In conclusion, we have proposed a traffic model based 
on packet routing strategy aiming at efficient paths in- 
stead of the shortest paths. This work may be useful 
for designing communication protocols for highly hetero- 
geneous networks since the present strategy can sharply 
enhance the throughput without any increase in its al- 
gorithmic complexity. The traffic rate process generated 
by this model exhibits critical self-similar scaling behav- 
ior with exponent ~ I, which implies the l//-type scaling 
in the power spectrum and the long-range correlated be- 
havior in a wide range of scales. The scaling behaviors of 
the present model is in good accordance with the empir- 
ical data, thus this model may have got hold of some key 
ingredients of the underlying mechanism of real traffic. 
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